function u = f0(x)
u=exp(-20.*(x-2).*(x-2))+exp(-(x-5).*(x-5));
end

